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Abstract: We compute the two- loop crossed six-line vertex master integral with two 
massive lines in dimensional regularisation, and give the result up to the finite part in 
D — 4. We apply the differential equation technique, and focus in particular on the purely 
analytical calculation of the boundary condition which we derive from a three-fold Mellin- 
Barnes representation. We also describe how the computation of the boundary condition 
is used to derive three non-trivial relations among harmonic poly logarithms of the sixth 
root of unity. 
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1. Introduction 

The computation of higher order perturbative corrections to heavy-to-light currents has 
recently been an active field of research. The analytical calculation of two-loop QCD 
corrections to differential semi-leptonic b — > u decays has been carried out simultaneously 
by several groups [1—4]. However, one of the occurring master integrals, namely the crossed 
six-line integral with two massive lines, has not been calculated purely analytically by any 
of the aforementioned groups. 

In this article we close this gap by reanalysing the master integral in question. We 
rederive the result through order 0((D— 4)°) by purely analytical steps, thereby confirming 
a suggestion in terms of transcendental constants for the boundary condition of the finite 
part in D — 4. During the course of the calculation we also found three non-trivial relations 
among harmonic polylogarithms (HPLs) of the sixth root of unity. 

This paper is organised as follows. In section [|| we define the kinematics and our nota- 
tion and summarize the final result. In section || we elaborate on the boundary condition 
which we derive from a three-fold Mellin-Barnes (MB) representation. In section || we show 
how to derive from the MB expressions three non-trivial relations between HPLs of the 
sixth root of unity. We conclude in section [|. 

2. Definitions and results 

In the following we consider the two-loop crossed six-line master integral depicted in Fig. [l|. 
We work in dimensional regularisation with D = 4 — 2e dimensions which is used to 
regularise both UV and IR divergences. Moreover, we define 



x = 



(pb - p) 2 + iy 



(2.1) 
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Figure 1: The two-loop crossed six-line master integral with two massive (double) lines of mass 
nib, and four massless (single) lines. The horizontal line has incoming momentum pb with p\ = ml, 
and the sloped line has outgoing momentum p with p 2 = 0. 



where +ir/ stems from the +irj prescription which we tacitly assume to be included in the 
propagators. Defining also the symbols 

our starting expression reads 
[dki] [ [dk 2 



(2.3) 



[{k 2 + Pb) 2 ~ ml] (k 2 + V? [(h + p b ) 2 - ml] (k 2 - h + pf 
1 

X (k 2 - h) 2 k\ ■ 

The integral in question was first given in [5] and subsequently used in [1-7]. Also the 
cases of six massless lines [8-13], as well as one [14], two [15,16], and four massive lines [17- 
19] are available in the literature. In the latter case the topology has two master integrals. 
Whereas the purely massless integral reveals a closed form in terms of hypergeometric 
functions of unit argument, the cases with massive lines are solved order by order in e by 
means of the differential equation technique [20-26]. 

The x-dependence of the integral / is obtained by taking the derivative of / w.r.t. x 
and plugging in the result of a Laporta reduction [27-29]. This yields a linear combination 
of the integral / itself and 16 other (fewer-line) master integrals, the latter are known 
to sufficiently high orders in e [5]. The solution of the differential equation yields the 
x-dependence except for the boundary condition at x = 0. 

Contrary to the case of four massive lines [17], the boundary condition at x = of 
the present integral cannot be inferred from the condition that the integral be regular at 
x = 0, but has to be calculated explicitly by means of other techniques. We summarize 
here our analytical result and postpone the derivation of the boundary condition at x = 
to section 0. Through order 0((D — 4)°) the result is 
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The value of the coefficient function at x = is one of the main new results of the 
present article. It confirms a conjecture given in [1, 3] which relies on the PSLQ [30] 
algorithm. However, we emphasize that all results in the present article were derived solely 
by analytical steps, and no fitting by means of PSLQ or other methods was involved. 

3. Computation of the boundary condition 

We now turn our attention to the computation of the boundary condition at x = 0. Al- 
though the coefficient functions cW are regular at x = 0, this property cannot be used 
in our case to pin down the integration constant from the differential equation method. 
Hence the boundary condition has to be computed by other means. Here we choose the 
Mellin-Barnes technique [31,32]. We start from the integrand in Eq. (|2.3j ), introduce five 
Feynman parameters to combine the propagators and subsequently integrate over the loop 
momenta. At this stage we set x = and carry out two of the Feynman parameter inte- 
grations. After applying analytic continuation formulas for hyper geometric functions [33] 
we arrive at 

I{x = 0) = -Si T 2 (l - e) T(l + 2e) (m 2 ) -2-2 * 
, 1 
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where u = 1 — u and similar for y and 2. Since each of the hypergeometric functions 
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possesses a one-dimensional Mellin-Barnes (MB) representation according to 

+ i oo 

P , , s T(c) fds T(- S )T(a + S )T(b + S ) s 
2 F l( a, b, c, -z) = J — z , (3.2) 

—i oo 

we can decompose each of the above terms via a three-dimensional MB representation and 
subsequently perform the integrations over u, y, and z in terms of T-functions. In one of 
the terms we can apply Barnes' second lemma. This yields 
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We preserve the order of the terms throughout this section. In the above equation ( |3.3; ) 
the contour integrals in the complex plane can be chosen along straight lines parallel to 
the imaginary axis, i.e. the real parts ki along the curves are constant. According to 
Refs. [32,34,35], these real parts, together with the parameter e, must be chosen in such 
a way as to have positive arguments in all occurring T-functions in order to separate left 



poles of T-functions from right ones. Therefore the four terms in Eq. (3.3) are regulated, 
respectively, for [36] 

5,5,1,9 . 

€ = -l28' kl = -8> k2 = ~V fe3 = -16' (3 ' 4) 
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e = , k\ = , ko = , fca = , (3.5) 

64 ' 32 ' 4 ' 16 ' V 1 

e = , k\ = — , ko = , (3.6) 

80 ' 5 40 ' v ' 

e = , fci = , ko = , k-i = . (3.7) 

768 384 ' 128 ' 192 v ; 

We perform the analytic continuation to e = with the package MB [36], which is also used 

for numerical cross-checks. We then apply Barnes' lemmas and the theorem of residues on 

the multiple Mellin-Barnes integrals, and insert integral representations of hyper geometric 

functions as well as ^-functions and Euler's l?-function where appropriate. The final result 

reads 
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The calculation of all but the second term in the above equation is more or less straight- 
forward. The second term, however, requires more effort. In this term, at the stage when 
all integrations have been carried out, we were still left with terms that contain HPLs of 
the sixth root of unity. Only after the application of the relations given in section || we 
were able to write the result in the form Q3.8]). 



4. Relations between HPLs of the sixth root of unity 

Higher transcendental functions that have as arguments powers of the sixth root of unity 
have been investigated at several places in the literature [37-42]. We introduce the notation 

X ± = l± Z -VS = e^s (4.1) 

from which follows immediately A+ + A_ = A + A_ = 1. The relations among harmonic 
poly logarithms of weight four that we found during the course of the calculation of the 
second term in Eq. (|3.3|) read 

- 5HPL({-3,1},A+) -23HPL({3,-1},A+) - 12 HPL({2, -1, -1}, A+) 

+6HPL({2, -1, 1}, A+) + 6HPL({2, 1, -1}, A+) 

_ 6 Li 4 (I)-Iln 4 2 + ^ 2 ln 2 2-^^ + C .c. =0, (4.2) 
K2J 4 38880 v ; 
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10HPL({-3, 1}, A+) + 73HPL({3, -1}, A+) - 18HPL({-2, -1, 1}, A+) 
-18HPL({-2, 1, -1}, A+) + 24HPL({2, -1, -1}, A+) + 21 HPL({-2, 1, 1}, A+) 

+21 Li 4 (k) + - In 4 2 - -vr 2 In 2 2 + + c. c. = , 

V2j 8 2 38880 
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- 5HPL({-3, 1}, A+) + 9HPL({3, -1}, A+) - 6HPL({-2, 1, -1}, A+) 
+3HPL({-2, 1, 1}, A+) + 4HPL({2, 1, -1}, A+) 

+ Li 4 (I) + lln 4 2-^ 2 ln 2 2 + I^ + C .c. =0. (4.4) 

Relations similar to the above ones were investigated systematically in [37]. We find it 
nevertheless instructive to give a few comments on their derivation. 

The first relation, Eq. ([4.2|), can be derived quite easily with the package HPL [43,44] 
by applying argument transformations to harmonic poly logarithms and exploiting the fact 
that A+ = 1/A_ and A+ = 1 — A_. 

The second one we obtained in the following way. In the second term in Eq. fl3.3j ), after 
the analytic continuation to e = 0, we cannot proceed until the end by merely applying 
Barnes' Lemmas and corollaries thereof. One way of performing the calculation in this 
case is to introduce new Feynman parameters as integral representations of higher tran- 
scendental functions (Euler's 5-function, ^-functions, hyper geometric functions) in order 
to be able to perform all MB integrations. Subsequently, one has to integrate over the 
newly introduced parameters. In our case, there are three such Feynman parameters to 
be introduced. Note that the three newly introduced Feynman parameters are different 
and independent from the three- fold integration in ( |3.lD which lead towards ([T^) 1 . At the 
stage where only one of these new Feynman parameters, say t, is left, we computed the 
occurring integrand twice by making two different choices of integrating by parts certain 
terms. The two results are then equated. 

The third relation, alike the second, also comes from the second term in (|3.3j), and also 
from the stage in where only one new Feynman parameter t is left in the integrand. Here 
we consider the four integrals 



t 

b 

j2 = ] it W-e)Mi zl±3 



o 



Ja = dt 



, LMirT+p]""'-'^ 



t 







1 We tried indeed to compute the second term directly from (3.1), without introducting MB integrations 



We succeeded to compute the simple pole, but not the finite part. 
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J 4 = dt V + J t . (4.5) 

o 

We compute each of them separately and subsequently equate them according to 

Li 2 (^ 1 _ t t + t2 ^j - 1) ln(l - 1 + t 2 ) + Li 2 (t - t 2 ) - Li 2 ^_ * +<2 ) = , < t < 1 . 

(4.6) 

Let us briefly outline the computation of one of the above integrals, namely J3. We first 
observe that ln(l — t + t 2 ) = ln(l — A+i) +ln(l — A_i). Therefore, we only need to compute 
the integral 

^(izhrt 2 ) 

dt V + / , (4.7) 



which we consider as a function of the complex variable A+. We first integrate by parts in 
order to be able to differentiate the polylogarithm. We then write 1 — t + t 2 = (1 — A+i)(l — 
t/\ + ). In the terms containing ln(l — t/\ + ) we subdivide the interval from . . . A + and 
A + . . . 1 due to the branch cut of the logarithm. We then take the derivative w.r.t. A + , 
which allows us to perform the integration over t. We finally integrate the result w.r.t. A + . 
The determination of the remaining additive constant is trivial. The computation of the 
other Jj proceeds along the same lines. 

Let us conclude this section by giving a few more relations between poly logarithms 
which turned out to be useful during the calculation. See also [40,45]. 

Li3a+ , +) ^ <3+ g M3)+j ( Sln2( 3 ) + ^-ftMi)), 



HPL({-2, 2}, A+) + HPL({-2, 2}, A_) 



1097T 4 vr 2 In 2 (3) $(-|,2,i) 2 



19440 



150 
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Jfc 
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$(-±3, ±) = Sv^Im 



50 

(4.9) 

(4.10) 
7rln(3) 10 /vrx 

(4.11) 
(4.12) 



All these relations also allow us to write the integrals J« in the following simple form 
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j 1397T 4 vr 2 ln 2 (3) 2<S>(-\,2,\) 2 6vr ln(3)$(-i , 2, |) 2vr <D(-i , 3, ±) 

2 1944 25 25 25\/3 5^/3 

(4.14) 

_ 989vr 4 7r 2 ln 2 (3) 2vr ln(3)$(-l, 2, |) 2vr $(-1, 3, |) 

3 9720 15 5^3 5V3 1 ' ' 
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(4.16) 

It is almost needless to say that also in this section all formulas were derived by purely 
analytical steps. 



5. Conclusion 



We presented the result of the two-loop crossed vertex master-integral with two massive 
lines. Our result represents the last missing piece in the fully analytic computation of the 
two-loop matching coefficients in heavy-to-light decays. The result was obtained by means 
of the differential equation technique. The boundary condition at x = cannot be inferred 
from the condition that the integral be regular at x = 0. Hence we computed it explicitly 
from a three-dimensional Merlin-Barnes expression. 

We emphasize again that all our results have been obtained by purely analytical steps, 
and no fitting of numerical values to transcendental constants by means of PSLQ [30] or 
other methods has been used. 
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